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Abstract: This paper gives a derivation for the large time asymptotics of the n-point density function 
of a system of coalescing Brownian motions on R. 


1. Introduction and statement of the main result 


The single species reaction-diffusion systems A + A ^ A (coalescence) and A + A ^ Q (annihilation) 
have been studied extensively in recent times [?^>| | 7||H| |n|ll()lllllll2llld| . A common interest in these papers 
is the departure from mean field behaviour and the computation of exact long-term asymptotics for the 
particle density. 

Recently, it was predicted in m that the large time asymptotics for the probability Pt{n,AV) of 
finding N particles in a fixed volume AV: 


Ptin,AV) 


t 2 4 

(if)”(lnt) 


d= 1 
d = 2 
d>2 


( 1 . 1 ) 


Note that the predictions for d > 2 agree with mean field behaviour. The second part of the exponent 
in d = 2 reflects multi-scaling, or deviation from linear scaling. In d = 2 the multi-scaling is manifested 
in the second logarithmic term. This type of scaling is indicative of particles being anti-correlated m 
These predictions were obtained by use of the dynamical perturbative renormalization group methods 
in a field theoretic setting |^. The setting here is for finite rate reactions, leading to annihilations of 
randomly walking particles on a fixed lattice. After renormalization the large time limit rate tends to an 
instantaneous reaction mu. Moreover the predictions carry over to the coalescing case, since they have 
the same effective field theory BIHl. 
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The aim of this paper is to verify the conjectures of na in dimension d = 1 , where detailed proba¬ 
bilistic tools are available. We consider a system of coalescing Brownian particles on the real line. Each 
particle evolves independently until it collides with another particle, at which time the two colliding par¬ 
ticles instantaneously coalesce into one. The n-point density function is defined, for distinct j/i,... ,2/„, 

by 

P [there exist particles in dyi ,..., dyn at time t] = Pnidh, ■ ■ ■ ,yn\t) dyi ... dyn- 

Note that pn depends on the initial particle distribution. The existence of such a density is dicussed in 
the appendix. 

Our main result can be stated precisely as the following asymptotic: 

n \y^-y 3 \ (1-2) 

l<z<j'<n 

where a{n) = n /2 -|- n{n — l)/ 4 . This asymptotic has the meaning that the left hand side is bounded 
above and below by constant multiples Ci, C2 of the right hand side. The upper bound is not actually an 
asymptotic, in that it holds simultaneously for all initial conditions, for alH > 0 and for all \yi\ < 
with a constant C2{L,n) < 00 depending only on n and L. The lower bound for all t > to and all 
\yi\ < Lt^!'^^ with a constant Ci(n, L, to) > 0 that depends on n,L,to and also on the initial condition. For 
the lower bound, the initial condition must satisfy a mild non-degeneracy condition; in particular it holds 
for deterministic initial conditions provided the gap between successive particles is bounded, but also if 
we assume the set of initial positions of particles {Xq : i > 1} is non-zero, translationally invariant and 
spatially ergodic, in the sense that the distribution of {Xq : * > 1} n (—00, 0] and {Xq : i > 1} n [T,oo) 
become independent as L —> 00. 

The key tool is the Karlin-McGregor formula for the non-coincidence probabilities for Brownian mo¬ 
tions. Useful upper and lower bounds on this transition density, which already display the key anomalous 
scaling term , are developed in section 0 by exploiting a representation known as the Harish- 

Chandra-Itzykson-Zuber formula (developed for random matrix problems). 

The empty interval method, and its generalizations, have been used to derive expressions for higher 
order correlation functions in [ 70 ^ for one-dimensional systems with instantaneous reactions. Large 
time asymptotics for the n-point correlation function density for the coalescing case with Poissonian 
initial conditions are found in [H], while the corresponding results for the annihilating case are given in 
[3. It was shown in this special case that the n-point density correlation function for the two systems are 
the same apart from the amplitude. This set of exact results was used to test the predictions mi in M 
The large time scaling of the formulae for the n-point density functions given in [HEl are not obvious as 
they involve a large combinatorial sum of terms with alternating signs. 

The problem of deriving rigourously the logarithmic corrections cu to the mean held answers in 
dimension d = 2 remains open. It would also be interesting to hnd out if there is a natural multi-fractal 
interpretation of the multi-scaling. Another simple system for which RG calculations predict multi-scaling 
in the stationary state is the system of aggregating massive point clusters with stationary source of light 
particles. This system is relevant to turbulence, see na It would be interesting to generalize the methods 
of present paper to prove multi-scaling for such cluster-cluster aggregation. 


2. Proof of the main result 


The proof is based on the following two lemmas. The hrst is a bound on the Karlin-McGregor formula 
for the transition density for non-intersecting Brownian motions 1^ . Fix x = {xi,... ,Xn) e R" and let 
(X“‘ : i = 1 ,..., n) be independent Brownian motions with Xg ‘ = Xi. For y £ R", denote by Gf^{x, y) 
the density of the probability measure 

P [Xf‘ £ dyi and the paths (Xf’ : s £ [ 0 , t]) are non-intersecting for i = 1 ,..., n.]. 
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Then the Karlin-McGregor formula is 


qKM 


{x,y) 


Gt{xi,yi) ... Gt{xi,yn) 


Gtixn^yi) ■ ■ ■ Gtixji^yn) 


( 2 . 1 ) 


where Gt{a,b) = is the one dimensional Brownian transition density. Note that the 

function Gf^{x,y) is the transition density for an n-dimensional Brownian motion killed on the set 
= yj}. One can check directly that the determinant expression ( 12.111 satisfies the heat equation 
with zero Dirichlet boundary conditions on this boundary, and that {x,y) = 5x=y 


The following lemma, proved in section 0 gives usable bounds on Gj^^{x,y). 


Lemma 1 For all xi <...< Xn, yi <■■■< yn and t > 0 

Gf^{x,y) 


, f f Gtixi^ yn—i-\-\) < 


i=l 




_i:2\ — Gt (Xj, J/i), 


( 2 . 2 ) 


i=l 


where and A{x) is the Vandermonde determinant defined by 

A{x) = A{xi,...,Xn) = (Xi-Xj). 


The second lemma, proved in section ^ gives a simple upper bound on the n-point correlation 
function, which reflects the intuition that particles should be anti-correlated, in that the presence of a 
particle in dx decreases the likelihood that there is a particle at another point dy. 


Lemma 2 For any initial distribution of particles, the n-point density function satisfies 

Pn{yi,---,ynG) < ( 7 rt)""G jgj. ally- g R and t > 0. 


The lower bound in the asymptotic follows quickly from the lower bound on the Karlin- 

McGregor formula (E3. Indeed, list the set of initial positions of particles as {X^ : f > 1} and let 
be the event that there exist particles , • ■ ■, X^" satisfying 

Xl'^ &[2kt^/‘^,(2k + l)t^/'^\ for yfc = l,...,n. (2.3) 

Then, applying the Karlin-McGregor formula between the points Xq^ , • • ■, N’q" and yi,..., 2 /„ we find, for 
t > 1, 


P [there exist particles at dyi ,. •., dyn at time t] 

> E [l(I2o(0) Gf ^ {{XIV. .., Xy), (2/1, ..., 2/„))] 


> Cr,E 


l{Q^{t))A{XlH-^l\ ..., Y[ Gt(X*^ 2/„_fc+i) 




A{yt dyi... dy„ 


>C{n,L)t ”GR[j7o(<)]Z\(2/t dyi... dy„ 


where we have used (lOl and \yj\ < in the final inequality. We have also used C{n, L,. . to denote 
a finite non-zero quantity, depending only on the quantities listed, but whose exact value is unimportant 
and may change from line to line. It remains only to bound P[i7o(t)] from below, independently of 
t > to- This clearly holds under the two sets of assumptions described in section ^ and in particular for 
Poissonian initial conditions. 
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For the upper bound in the asymptotic we estimate the probability that there exist particles 
at dj/i,..., dyn at time 2t by conditioning on the set {X^ : i £ N} of positions of the particles at time t. 
For the desired particles to exist at time 2t, one of the events 

= {Xll G dyk and the paths : s £ [t, 2t]) do not intersect for /c = 1 ,..., n} , 

for some < ... < must occur. Applying the Markov property at time t and the upper bounds in 
Lemma n we find, for all t > 0, 

P [there exist particles at dyi, ■ ■ ■, dyn at time 2t] 


< Cn E 


A{yt dyi... dy„ 


A{X\H-P\...,Xl-t-P’^) R Gt{X\\yk) 

< C{n) Gt{x, y)A{xt~'^/'^)pn{x] t)dx'^ A{yt~^/'^) dyi... dyn see (IB.Ill 

< C{n, L) A{yt~^/‘^) dyi... dyn- 


The final inequality uses the substitution x 


xt}!"^ and the bound from Lemma [3 


3. Proof of Lemma m the Karlin-Macgregor bounds 


First note that Brownian scaling implies that Gf^{x, y) = t Gi’^ {xt yt and we need only 

prove the lemma for t = 1. Factoring out common terms from the rows in EU we find that 

n 

Gf^(x,y) = (27r)-"/2det(F;(x,y)) 

i=l 

where E(x,y) is the matrix with entries exp(a;i?/j). We now use the Harish-Chandra-Itzykson-Zuber 
formula |21, which is widely used in random matrix theory. This states that 

det(£'(a;,y)) = c„Z\(a;)Z\(y) [ exp {tr{UXU'^Y)) p{dU), 

Ju{n) 

where tr(A) is the trace of a matrix A, fi{dU) is (normalized) Haar measure on the unitary group U(n) 
and X, Y are the diagonal matrices with entries xi,... ,Xn and j/i,..., A short proof of this formula 
is found in ^ (appendix A.5). 

The fact that U is unitary ensures that the function E{U) = tr([/AT7^y) is real valued (recall that 

W is the conjugate transpose of U so that UW = U^U = I). By the compactness of U(n), F{U) achieves 

its maximum and minimum values and we claim that 

n n 

min E{U) =y^ Xn-i+iVi, max E{U) =y^x^yi. (3.1) 

UGUM t/eU(n) 

2=1 2=1 

Bounding the integral in the Harish-Chandra-Itzykson-Zuber formula using these maxima and minima 
in the integrand leads to the bounds 12 . 211 . 

It remains only to prove the claim EU, for which we will find all the stationary points of the smooth 
function E on the manifold U(n). A matrix V is in the tangent space to U(n) at the point U precisely 
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when UV^ + VW = = 0. The first derivative of F in the direction of V in the tangent space 

is given by 

DvF{U) = tr{UXV^Y) +tr{VXU'<Y) 

= ti{V^YUX) +\:T:{yXU^Y) (since tr{AB) = tr{BA) for Hermitian B) 

= tr(C/yty[/xC/t) + tr(yXC/1'y) (since tr(C/^C/t) = tr(yl)) 

= -tv{VU'^YUXU^) + tviVU'^UXU^Y) (using = -VU^) 

= tr {VU^ [UXU^,Y]) 

(where \A, B\ = AB—BA is the commutator). The map V VU^ bijectively maps the tangent space onto 
the Lie group Lie(U(n)) defined by {IT : + W = 0}. It is straightforward to check that tr(lTA) = 0 for 

all W £ Lie{XJ{nj) implies that ^ = 0. The vanishing of DvF{U) therefore implies that \U XU\Y\ = 0. 
This in turn implies, since Y is diagonal with distinct eigenvalues, that UXW is a diagonal matrix. The 
matrix UXU^ has the same distinct eigenvalues as X, so that the entries on the diagonal oi UXW must 
be a permutation of xi,... ,x„. Hence the value of F{U) at a stationary point is for some 

permutation tt. We must show that this value is maximized when tt = Id the identity permutation and 
minimized at the involution 7 r(i) = iV + 1 — L Argue by contradiction and suppose that tt — > ^Tr(i)yi 
is maximized over permutations at some ttq ^ Id. Then there exist i < j for which 7 ro(*) > 7 ro(j). Let 
be the transposition permutation, swapping the jth and fcth elements. Then an increased value is 
found by taking the composition tti = o ttq permutation, since 

n n 

^ ^ ~ ^ ^ ^'Ko(i)yi (^7ro(j) ^'^o{j)')iy3 ~ V'i-) ^ 

i—\ i—1 

using the ordering xi < .. . < Xn^ yi < • ■ ■ < Un- The minimum follows from a similar argument. 


4. Proof of Lemma n the n-point density bound 


We start with a simple construction of coalescing particles as follows. Fix xi < X 2 < • ■ Take I.I.D. 
Brownian motions (L* : i > 1) and set Xf* = Xi + Bl, for t > 0, the non-interacting Brownian paths. 
Then define Xy = Xy and inductively for i >2 

X^i ^ f for t < T(xi,x^-i), 

‘ \X^'~^ for t > r(xi,a;j_i), 

where T(xi,Xi-i) = infjf : Xf* = X^'~^} = inf{t : Xy = X^'~^}. The path t Xf gives the position 
of the particle that started at x. This construction shows that the process starting from infinitely many 
particles can be approximated by the process starting with finitely many particles in an increasing way. 
It will therefore be sufficient to bound pn over all finite starting positions. We write Pxi,x 2 ,... fo indicate 
the starting positions. 

The fact that pi{y) < ( 7 rt)“^/^, for any initial particle configuration, is well known. In section 9 of 
m it states that the time reversal duality formula for coalescing Brownian motions for a < b is given by 

Pxi,...,xm [{^r ■i<m}n[a,b]^9] = Pa,b [[X“,X{'] n {xi,.. .,Xm} ^ 0] . 

We can take the supremum of the above duality over all starting configurations to obtain 

sup P,,[{Xy : z > 1} n [a,5] ^ 0] = Pa,b [X“ < Xf] 

Xi 

which can be explicitly calculated and is bounded by ( 7 rt )“^/^(6 — a). We give an elementary approach 
to this bound that avoids duality in Appendix 0 
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For the extension to pn given in Lemma |21 we fix oi < 5i < 02 < ... < a„ < Define 

= {there exists a particle in [oj, bj\ at time t} . 

We shall show by induction that 

n ... n (4.1) 

i—1 

Decompose f2i as the disjoint union f2i = Ukf^i.k where 

^i,k = {^t' < for * = 1 ,..., fc — 1 and G [oi, 61 ]} . 

Then conditional on ^ the processes : i > k + 1) evolve as coalescing Brownian motions with an 
extra lower absorbing boundary along the path t X^*‘. The construction above shows that this extra 
absorbing path only lowers the probability that any of the paths {X^' : z > fc + 1 ) reach the intervals 
[ 02 , 62 ],..., [a„, 6 „]. Thus, applying Bayes formula, 

m 

Px\,...,Xm n . . . n f2ri] — ^ ^ Pxi,...,Xrri [^2 H . . . H f2ri I ^l^k] Pxi,...,Xm. [‘^ 1 ,^] 

k^l 

m 

^ n ... n f2n] Pxu...,Xm [^hk] 

k^l 

m 

< [f ?2 n... n Qn] Pxu...,x^ [f 2 l,fc] 

k^l 

= Pxu...,x^ [122 n... n 12 „] [ 12 l] . 

and the result 631, and hence Lemma El follows. 


A. An elementary approach to the one point density 


This elementary approach may be of some interest for systems that do not have a duality relation. 
Fix a € [0,1] and let p(xi,... ,Xn) be the probability that, starting some coalescing Brownian motions 
xi,... ,Xn, there is a particle at time t in the interval [0, a]. We aim to show, by induction on n, that 

p(xi,... ,x„) < (A.l) 

for some constant C < 00 . We suppose < ... < Xn+i- Condition on the path t —> of the path 

started at Xn+i- We can consider three cases: if < 0 then we are sunk since no other path can 

ever enter [0, a] due to coalescence; if G [0, a] then we are done; finally if > a then we still 

need some of the particles from xi,..., to enter [0, a], but there is an absorbing boundary along the 
path of This absorbing boundary can only lower the chance of getting particles where we want. 

This splitting of possibilities leads to 

p(xi,..., Xn+i) < C't-^/2aP[A("+i) > a] + P[A(”+i) G [0, a]]. (A.2) 


It is straightforward to find cq < cxd so that for a Brownian motion (A*), and for L > 1, t > 0, 


PjXt G [0, g] I Aq = x] 
P[Xt < a I Ao = x] 


< coLat 


whenever jxj < Lt 


Suppose now that Xi G [—for all i. Using this bound the induction argument works with the 
choice C = cqL. In particular we have shown that 


sup p(xi,..., Xm) < coLt 

Xi|<Ltl/2 
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Applying this with the choice L = 0(-\/log(t)), and using a simple bound that any particle starting 
outside [—can reach the intervals [ 0 ,a] (we need only control particles starting exactly at 
by coalescence) we find that 

supp(a;i,.. .,Xm) < C't“^/^(log(t))^/^a 

Xi 


where the supremum is over all possible initial configurations. One can remove the unwanted logarithm 
term by a blocking argument, but this is already longer than the elegant duality argument, which has 
the extra advantage of achieving the optimal constant. 


B. A note on the existence of pn 


We briefly consider for which measure t) acts as a density. List the (disjoint) positions of the particles 
at time t as {XI : i > 1}. Define a random measure on R" by 

E W,...,Ay) 

where S{x) is a point mass at x. Let Cn{S) denote the cube [0, i^)”. The cubes k + C'„(2“"*), for k g 
(2“™Z)", partition R”. Take A C R" open. For an initial condition with only finitely many particles, 
the measure /i" has only finitely many atoms which do not lie on Ui^j{xi = Xj}. Therefore we have the 
increasing limit 

hm ^ lipnk + Cni2-^))>0)=p^iA), 
fee{2-^Z)^ 

where the sum over k is restricted to those terms for which fci,..., are distinct. The event /r"(fc + 
C'„(2“™)) > 0 can be rewritten as [kj, kj + 2“™] n : i > 1} 0 for j = 1,..., n. Taking expectations 

and applying 14.111 we find that E[fdp{{] < (27rf)“”/^ J fdx when / = 1(A). By an approximation 
argument this also holds for all measurable / and by monotonicity it holds for any initial distribution of 
particles. We may then define p„(a;; t) = p{xi, ..., a;„; t) as the density of the measure E[p]{], that is the 
Radon Nicodym derivative with respect to Lebesgue measure, so that 



f{x)pn{x]t)dx = E 


H /(F‘. At) 


(B.l) 


The Lebesgue differentiation theorem implies that p(x; t) = lim^^o S "A[^”(x + C„((5)] for almost all x. 
Choosing xi,... ,Xn disjoint we obtain 

p{xi, ...,Xn;t) = lim S~'^E [Nt{[xi,xi + (5])... iV* ([a;„, + S])] 

o—*0 


where Nt{A) is the number of particles inside A at time t. However we claim we may replace A(([a;i, a;i+(5]) 
by the indicator I([a;i, Xi + 5] contains a particle) in this limit. The error when doing this occurs if there 
are two or more particles in [xi,xi +5] and is dominated by E[p^'^^{{xi,xi,X 2 ,..., x„) + Cn+i (<^))]- This 
error is therefore of order (5"+^, so that replacing Nt ( [xi,Xi + 5]) in this way for each z = 1 ,..., n we reach 
the usual definition of Pn{x\ t) as, for distinct xi,... ,Xn, 

p{xi^... ,Xn]t) = lim 6~'^P [all the intervals [a ;^,Xi + <5] for z = 1 ,...,n are non-empty]. 

( 5^0 
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